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I.  INTRODUCTION 


-.  /  /  </  i 

V  ■  /  f-  .«■  o 

^  '  ■'£  f  V 

\  v 


•  In  the  past  several  years,  there  has  been  widespread  interests 
in  a  variety  of  wave  generation  devices  based  on  un-neutral i zed  beams 
of  relativistic  or  mildly  relativistic  electrons.  Such  devices  are 
potentially  capable  of  producing  high  power,  with  high  efficiency,  in 
the  centimeter,  millimeter,  and  sub-millimeter  SDectral  ranges.  An  , 


example  is  the  electron  cyclotron  maser,  which  has  produced  AlOO  kW 
at  centimeter  wavelengths,  with  efficiencies  ranging  from  30%  -  50%. 

The  gyrotron  mechanism  involves  azimuthal  phase  bunching  of  the 

electron  beam,  which  occurs  during  resonant  interaction  of  the  Doppler- 

_ 

shifted  beam  cyclotron  mode  with  a  waveguide  or  cavity  modejZJ-O^'rr 
The  bandwidth  of  the  gyrotron  mechanism  is  typically  of  order  1-2%.  , 


For  various  military  and  communications  applications,  it  could 

be  of  great  interest  to  develop  free-electron  amplifiers  with  wide 

operating  bandwidths,  AwJLuiJ'  In  this  report,  we  present  a  linear 

analysis  of  a  new  principle  for  wideband. ampl ification.  The  basic 

configuration  is  the  use  of  an  annular  electron  beam  drifting  in  a 

JL^ 

dielectrically-loaded  waveguide  In  such  a  waveguide,  the 

normal  guide  modes  have  phase  velocities  asymptotic  to j  c/^e7 where 
c  islthe  speed  of  light  and — tU— vs  the<di electric  constant^ (assumed 
to  be  isotropic)  of  the  waveguide  lining. "'-The  system  is  assumed  to  be 
immersed  in  an  axial  magnetic  field*  sufficiently  strong  to  render 
the  motion  of  the  beam  electrons  one-dimensional.  In  this  case,  — * — — 
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the  amplification  mechanism  depends  on  axial  charge  bunching  in  the 


Cerenkov  interaction  of  a  longitudinal  space-charge  wave  of  the  beam 
with  the  !yz  -  componant  of  a  TM  guide  mode^ 

The  interaction  we  envision  is  schematically  represented  in 


the  dispersion  diagrams  of  Fig.  3,  where 


p 

u>+  =  kV  ±  - — 

- 

V  £Y 

are  the  positive  and  negative  energy  beam  modes  and 


are  the  waveguide  modes  for  the  case  of  a  constant  isotropic  dielectric 
material  filling  the  waveguide.  As  can  be  seen  in  Fib.  3a,  where  s=l , 
the  beam  and  guide  modes  do  not  cross  and  therefore  four  distinct  pro¬ 
pagating  modes  are  expected  on  the  entire  k  spectrum,  with  their  dispersion 
characteristics  given  by  slight  modification  of  Eqs  (1)  and  (2)  due  to 
coupling.  On  the  other  hand  however,  Fig.  3b,  when  e  +  1  and  in  particular 
when  e>  (c/V)  the  dispersion  curves  do  cross  and  therefore  growing 
modes  are  expected  due  to  the  interaction.  Most  interesting  appears  to 
be  the  case  when  e  slightly  exceeds  the  value  (c/V)  ,  since  then  strong 
interaction  is  expected  over  a  broad  band  of  frequencies. 

The  above  discussion  applies  also  to  the  topologically  equivalent 
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case  of  a  waveguide  with  partial  dielectric  filling,  for  which  the 
dispersion  curve  is  represented  by  the  broken  line  in  Fig.  3b.  Since 

o 

this  configuration  is  more  compatible  with  e  ~  (c/V)  ,  the  analysis 
of  this  report  is  based  on  the  realistic  geometry  of  Fig.  2. 

In  Section  II,  we  present  the  general  linear  analysis  of  this 
system,  culminating  in  the  dispersion  relation  of  Eq.(ll).  In  Section  III 
the  dispersion  relation  is  solved  analytically  in  a  simple  limiting  case. 
In  Section  IV,  we  present  the  conclusions  of  this  analysis  and  suggestions 
for  future  work  on  this  concept. 
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II.  ANALYSIS 


A.  Physical  Model 

The  configuration  and  geometrical  parameters  are  defined  in 
Figure  2.  The  waveguide  is  a  perfectly  conducting  sylinder  of  radius 
r^.  The  guide  is  lined  with  a  lossless  dielectric  (of  isotropic 
dielectric  constant  t)  in  the  space  r^  <_  r  <.  r^.  The  electron  beam 
is  annular,  filling  the  region  r-j  <  r  ^  r~.  The  regions  0  $  r  <  r-j  , 
r2  <  r  <  r3  are  in  vacuum-  An  applied  axial  magnetic  field  Bq  =  BQez 
pervades  the  waveguide.  The  un-neutral i zed  electron  beam  has  purely 
axial  motion;  this  corresponds  to  the  limit  Bq  ■+•  °°,  so  that  transverse 
motion  is  suppressed  and  no  cyclotron  interaction  may  occur.  We  assume 
a  cold  uniform  annular  beam,  the  equilibrium  distribution  function  of 
which  is  given  by 

f0  (r»P2)  =  nQH  (r-r-j )  H  (r2-r)  5  (pz  -  YQmV )  (3) 

2  -1 

where  =  ( 1 -V  )  and  H(x)  is  the  Heaviside  step  function,  and 

2  2  r 

where  ojr  =  4irn  e  ly  m. 
b  o  J  o 

As  discussed  in  the  introduction,  when  V>  c//~  so  that  the  un¬ 
perturbed  beam  and  guide  modes  intersect,  we  expect  instability  (Fig.  4). 
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We  seek  to  derive  the  general  dispersion  relation  for  this 
instability,  in  order  to  obtain  quantitative  results  for  the  parameters 
of  interest  for  the  design  and  operation  of  a  device  to  act  as  a  broad¬ 
band  amplifier.  Such  parameters .indicated  in  Fig.  4,  are: 

•  The  optimal  frequency,  cu  ,  for  which  the  amplification 
rate  is  largest. 

•  The  corresponding  spatial  growth  rate,  6k j. 

•  The  cutoff  frequency,  a  ,  corresponding  to  the  modified 

waveguide  curve.  J 

•  The  upper  cutoff  frequency,  ou, ,  of  the  unstable  range 

•  The  choice  of  values  for  e,  V,  n  etc  for  which 

oj,-w  and/or  5k,  maximize. 

I  m  I 

To  derive  the  dispersion  relation  for  the  coupled  beam-guide 
system  of  Figure  2,  we  consider  TM  .modes  (i.e.,  no  azimuthal  de- 

U  U 

pendence) . 

The  wave  fields  are  given  by 

t  (r,z,0,t )  =  [ez(r)  ez  +  Ep( r)  er  ]eikz‘iut  ,  (4) 

3  (r,z,0,t)  =  BQ(r)  e9  e  lkz-iwt  ,  (5) 

where  E  (r)  is  the  solution  of  the  wave  equation 

h2(w,k,r)  | ( r )  =  0,  (6) 

with  E„  and  Ba  given  by: 
r  i 

1 1-2 


_d_ 

dr 


k  -e(r)  -SSj 
c 
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Er  (r) 


ik _ 

e(r)  aj2/c2-k2 


(7) 


! 


H 


0 


gj 

Ck 


(8) 


and  =  eEr,  HQ  =  Bg/y  ,  with  y  being  the  magnetic  permittivity. 
2 

The  quantity  h  (u),k,r)  is  the  plasma  dielectric  function 


2 

h  (w,u,r)  =  1  - 


2 

4tt  g 
m  c(r) 


dPz  fo(v’Pz) 
Y3  (eg  -  kvz)2 


(9) 


where  the  equilibrium  distribution  function  f0(r,pz)  given  by  Eq.3 
will  be  assumed. 

Using  Eqs.  (4)  -  (9),  the  dispersion  relation  is  obtained  by 
applying  the  boundary  conditions: 


Ez  (r4)  =  0  , 

[Ez(r)]  =  [B0(r)]  =  0  (r  =  rl’r2’r3)  • 

where  £  r )J  indicates  the  jump  in  quantity  y (r)  across  the  interface 
at  r.  The  seven  conditions  (10)  yield  the  dispersion  relation. 
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02a) 
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02b) 

02c) 


the  last  expression  being  the  dielectric  for  a  cold  beam,  f0~^(Pz  *  PQ) 

and  where  J  ,  Y  (I  ,K  )  are  the  ordinary  (modified)  Bene!  functions  of 
nr  m  nr  m 

order  m.  Equation  (11)  is  the  central  result  of  the  analysis.  It  must 
be  solved  numerically,  giving  the  wavenumber  (Re  k)  and  spatial  growth 
rate  (  -  Im  k)  for  real  iu. 
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III.  APPROXIMATE  SOLUTION  OF  THE  DISPERSION  RELATION 

An  approximate  analytic  solution  of  Eq.  (11)  can  be  obtained  in 
the  limit  where  we  consider  a  thin  annular  beam  at  the  dielectric  inter¬ 
face;  i.e., 

r^  -  r2  -*■  0;  r ^  =  r^r^  «  r^  .  (13) 

The  introduction  of  this  limit  is  done  not  just  for  the  sake  of 
the  simplification  and  tractability  it  introduces  to  the  equations,  but 
also  because  it  corresponds  to  that  of  the  most  efficient  operation.  In 
fact,  since  the  waveguide  modes  outside  the  dielectric  material  are  pro¬ 
portional  to  the  modified  Benel  function,  their  amplitude  is  maximum  near 
the  dielectric  interface,  hence,  it  is  desirable  to  place  the  beam  the 
closest  possible  to  the  dielectric  to  maximize  the  coupling  coefficient. 

Under  these  assumptions,  Eq.  (11)  can  be  simplified  to  yield  the 
dispersion  relation 


1 


Wz’ 

W? 


~VkrrV  -  VkEr2^  'Vk£r4i 


*  rZlh 


0  . 


(14) 


From  Eq.  (14),  we  may  extract  several  limiting  cases: 
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(a)  If  either  k ^  »  1  or  k^r^  =  kE(r4  -  r2)  «  1 ,  we  have 


J1(k£r2)  Y0  (k£r4)  -  Y1(kEr2)  JQ(kEr4) 
JAr2>  Yo  (kEr4}  -  Yo(kEr2>  J0(kEr45 


ctn  (k£r42) 


(15) 


gi ving 


1_  Il(ltor2^ 
ko  W? 


1^  ctn  (^r^)  +  r2i^2  "  ® 


(16) 


(b)  If  kQr2  >>  1 ,  then 


Il(kor2) 

W? 


2  kor2 


giving 


(17) 


r  •  b  ctn  (kEr42> 

0  E 


+  r21h2  =  0 


(18) 


Equation  (17)  is  not  valid  near  the  cutoff  coQ.  The  result  (18)  is 
indentical  to  that  obtained  in  planar  geometry;  this  is  the  limiting 
case  of  large  aspect  ratio. 


(c)  Near  the  beam  modes,  kEr42  <.t/2,  and  we  have 


TT 

2 


42 


r21h' 


The  cutoff  frequency  ajQ  (k  =  0)  is  obtained  from 


(19) 
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*1  •]  ( ^  ^ 
Jo(ajor2/c) 


(20) 


For 


eV 


42 


»  rr 


implying  ui  rg/'c  <  1 ,  we  find 


rrc 


2  z*r 


42 


TC 


For  £'2r42  ~  r2 


(wQr2/c  >,  1 ) ,  we  have 


(•21) 


jjq  =  2.405  c/r4  . 


(22) 


Near  the  beam  mode,  we  have  from  Eq.  (19)  the  form 


(w  -  kc3) 


1  TT  £ 


Lko (^’k )  '  2  kE^,k) 


jy +  er42 

m 


.  V 
21  V 


where  6  =  VQ/c. 


Taking  k  =  uj/Bc  +  5k,  where  |5k[  <<  ikU 


we  find 


(<5k)2Dg(oj,k)  =  , 


where 

Og(u,k)  «  ‘  +  £r42 


(23) 


(24) 


(25) 
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expanding  Dg  about  k  =  w/6c  ,  we  have 


(<$k): 


3D_ 

°g(u)’  tt]  +  3k1  (a)’ 


r2i 

Y3S2c 


2  » 


(26) 


which  can  be  solved  for  Sk(oj) .  The  maximum  growth  rate  is  obtained  at 

U)  =  ium  for  DJor,  w/8.)  =  0,  giving 
m  g  m  m  c 


m 


rr 

2 


3  1 

/eS2  -"l  - 


*  v^T 


(27) 


as  the  frequency  at  which  the  growth  rate  maximizes;  the  corresponding 
growth  rate  is  given  by  the  imaginary  part  of  6km,  where, 


r->iJj 

6k 3 

m  r2  c2 
r42c 


=  2rp  1 
2 


eB2Y3 


(28) 


From  Eqs.  (27)  and  (28)  we  deduce  an  important  constraint,  viz ,  that  the 
spatial  growth  rate  and  frequency  bandwidth  are  universely  related.  Large 
growth  rates  can  be  obtained  at  the  expense  of  limited  bandwidth,  and 
vice  versa.  To  see  this,  we  note  that  j(5km)3|  maximize  for 

e  >>  1 ;  s32  -*■  2  ,  (29) 


where  we  set  y  —  1  (i.e.,  the  beam  is  only  weakly  relativistic).  In  this 
limit,  we  find 


nc  3 
r42 


(30) 
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and 


?  ('  - 


(l  -  i  /t)  i 

r2l“b  TT  fe82  -  l)a 

A 

\  /  Y 

/42c!  c6! 

1/3 


(31) 


Thus  in  the  limit  e2  r ^  >  r2  (f°r  £  >>  1)»  we  have 


^  =  4b2 
<"o  V  £6L  - 


-  /T 


(32) 


under  the  conditions  (29)  for  which  J 5k  I  is  maximum. 

m 

On  the  other  hand,  the  bandwidth  6oj  is  given  by  the  upper  cutoff  fre¬ 
quency  co-j ,  where  Im<5k  =  0,  as 


,  .  .3 

-  u,  -  U)  „  - 

1  m  2y 


r42  e 


1/3 


(33) 


Thus  from  (28)  and  (33)  we  find  the  very  useful  relation 


<5u) 


2.18 


clm5km  ,u  e62  -  1 


(34) 


which  yields  a  maximum  for  5w  at 
£  ♦  2,  c82  ■*  1 ,  y  -*  /2  ; 


(35) 


the  regime  (35)  is  seen  to  be  the  opposite  limit  to  that  of  relation 
(29). 
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IV.  DISCUSSION  AND  CONCLUSIONS 


We  shall  briefly  discuss  the  consequences  of  relaxing  a  simplifying 
assumption  we  have  used  in  the  analysis,  viz_  that  the  beam  is  taken  to  be 
cold  and  uniform  in  density.  We  have  neglected  effects  of  temperature  and 
energy  shear.  Across  the  beam,  the  wave  field  is  given  approximately 
(for  r21  «  rx)  by 


SrT-  ko  h2  ^ 


Ez(r)  =  0  , 


where 
h2(r)  *  1 


iSSi  ffo(r’p)  dp 

m  J 


(36) 


(37) 


Y3 (w  -  kv) 3 

If  h(r)  is  a  smooth  function,  Eq.  (36)  has  the  WkB  solution 


E-(r)  *  h  exp  (k  /  h(r)  dr)  . 


./* 


(38) 


'We  may  estimate  the  effects  of  energy  spread  and  self  fields  (manifested  in 
the  potential  s(r),  by  taking  an  equilibrium  distribution  of  the  form 


/ fur '  e<>(r)  exP 

_  (p2/2m  -  e<p(r)  -  EQ) 2 

L  2  E|h  J 

(39) 


where  EQ  is  the  nominal  beam  energy  and  E^  is  a  measure  of  thermal  energy 
spread.  Then  the  quantity  w£/y3(w  -  kv)2  used  for  the  analysis  of  a  cold 
beam  is  to  be  replaced  by 
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where  quantities  subscripted  "zero"  are  evaluated  at  the  center  of  the 
beam  and  where 

e<5<|>  =  mo)2b0  r^/8  (41 ) 

is  the  potential  energy  drop  between  the  edge  and  the  center  of  the  beam. 
The  energy  shear  associated  with  the  self- fields  has  two  effects: 

(i)  The  effective  beam  density  is  reduced: 


(ii)  The  effective  energy  spread  is  increased: 
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Summary 

In  the  analysis  of  Sections  II  and  III  we  have  demonstrated  the 
validity  of  the  dielectric-loaded  Cerenkov  wide-band  amplifier  concept. 

We  have  obtained  approximate  relations  for  the  growth  rate  and  bandwidth, 
revealing  an  important  tradeoff  between  them,  cf.  Eq.  (34),  which  must  be 
considered  in  design  applications. 

An  important  extension  of  this  work  would  be  to  consider  a  nonlinear 
calculation  in  order  to  elucidate  the  saturation  mechanism  (which  we 
expect  to  be  beam  trapping)  and  ultimate  power  levels. 
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FIGURE  CAPTIONS 


Figure  1  -  Dispersion  curves  for  the  waveguide  mode  and  the  electron 
cyclotron  mode  in  the  electron  cyclotron  maser. 

Figure  2  -  Geometry  of  the  Broadband  Amplifier. 

Figure  3a-  Dispersion  curves  for  the  waveguide  mode  and  the  beam  modes 
when  z  -  1 . 

Figure  3b  -  Dispersion  curves  for  the  waveguide  mode  and  the  beam  modes 
when  the  waveguide  is  filled  with  a  dielectric  material  with 
c  >  (c/V)2.  The  broken  curve  refers  to  a  partially  filled 
waveguide. 

Figure  4  -  Wavenumber,  Re(k),  and  spatial  growth  rate,  Im(k)  for  the 
instability  (schematically).  The  solid  curves  indicate  the 
unamplified  propagating  modes. 
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